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Abstract 

The two-parameter Birnbaum-Saunders distribution has been used succesfully to model fa- 
tigue failure times. Although censoring is typical in reliability and survival studies, little work has 
been published on the analysis of censored data for this distribution. In this paper, we address the 
issue of performing testing inference on the two parameters of the Birnbaum-Saunders distribu- 
tion under type-II right censored samples. The likelihood ratio statistic and a recently proposed 
statistic, the gradient statistic, provide a convenient framework for statistical inference in such a 
case, since they do not require to obtain, estimate or invert an information matrix, which is an 
advantage in problems involving censored data. An extensive Monte Carlo simulation study is 
carried out in order to investigate and compare the finite sample performance of the likelihood 
ratio and the gradient tests. Our numerical results show evidence that the gradient test should be 
preferred. Three empirical applications are presented. 

Key words: Birnbaum-Saunders distribution; Censored data; Fatigue life distribution; Lifetime 
data; Likelihood ratio test; Maximum likelihood estimation; Monte Carlo simulations; Gradient 
test. 

1 Introduction 



Birnbaum and Saunders 



(|1969aUbr) proposed a family of two-parameter distributions to model failure 
time due to fatigue under cyclic loading and the assumption that failure follows from the development 
and growth of a dominant crack. This distribution is known as the two-parameter Birnbaum-Saunders 
(BS) distribution or as the fatigue life distribution. The BS distribution is an attractive alternative to 
the Weibull, gamma, and log-normal models, since its derivati on considers the b asic characteristics of 



the fatigue process. A more general derivation was provided by Desmond (1985) b ased on a biological 



mode 



and relaxing several of the assumptions made by 



Birnbaum and Saunders 



(1962ab. 



Desmond 



(|1986|) investigated the relationship between the BS distribution and the inverse Gaussian distribu- 
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tion and demonstrated that the BS distribution is an e qual- weight mixture of an i nverse Gaussian 
distribution and its complementary reciprocal; see also lBhattacharyya and FriesI (|1982l) . 

The random variable T is said to have a BS distribution with parameters a, (3 > 0, say BS(a, (3), 
if its cumulative distribution function (cdf) is given by 



F(t) = $0) 



t > 0, 



(1) 



where $(•) is the standard normal distribution function, v = p(t//3)/a, p(z) = z x l 2 — z x l 2 and 
a and (3 are shape and scale parameters, respectively. Also, j3 is the median of the distribution: 



F(p) = $(0) = 1/2. For any constant k > 0, it follows that kT ~ BS(a, krj] 
the reciprocal property holds for the BS distribution: T _1 ~ BS(a, see 
probability density and hazard ratio functions corresponding to © are given by 



. It is noteworth y that 



Saunders 



(1974). The 



f{t) = K{a^)t-*/ 2 {t + (3) 



exp 



2a 



t > 0, 



(2) 



and 



r(t) 



n(a, /3)t" 3/2 (£ + 0) exp{-r(t//3)/(2a 2 )} 
1 - $(v) ' 



t > 0, 



respectively, where K(a,ft) = exp(a 2 )/(2ay/2ir/3) and r(z) = z + z . The expected value, 
variance, skewness and kurtosis are, respectively, E(T) = (3(1 + a 2 /2), var(T) = (a(3) 2 (l + 5a 2 /4), 
p 3 = 16a 2 (lla 2 + 6)/(5a 2 + 4) 3 and /i 4 = 3 + 6a 2 f 93a 2 + 41)/(5a 2 + 4) 2 . A general expression 



Rieck 



19991) and is given by 



to obtain the moments of © was derived by 

' K p+1/2 (a- 2 ) + K p _ 1/2 (a- 2 ) 



E(T P ) = ^ 



2K 



1/2 



or 



where K v (-) denotes the modified Bessel function of the third kind and order v. As noted before, 
if T ~ BS(a,f3), then T' 1 ~ BS(a, /T 1 ). It then follows that E(T _1 ) = ^^(l + a 2 /2) and 
var(T _1 ) = a 2 [3~ 2 (l + 5a 2 /4). The BS distribution is positively skewed and the asymmetry of the 
distribution decreases with a. Figure \T\ plots the density and hazard ratio function for some values 
of a with (3 = 1. For applications of the BS distribution in r eliability st u dies s ee, f or example 



B alakri shnanetal 



Leiva et al. 



(2008D and 



Leiva et al 



(|2007|) . and for its use outside this field see 

12009D. 

The BS distribution has received significant attention over the last few years, 
discussed the shape of the hazard fu nction of the BS dis tr ibution. Improved freq uenti st inference for 
the BS distribution is discussed in 



Kunduetal 



(2008) 



Wu and Wongl (bo041) . Icvsneiros et all (120081) an d 



Guiraud et al 



Castil 



(I2009h and 
conditional duration model and 
based on the truncated life tests. 



o et al 



Lemonte et al. 



(|2007l |2008|) . while estimation using referenc e prior is presented in IXu and T ang (20 10|). Exten - 
sions of the BS dis tribu tion are proposed in Dfaz-Garcfa and Leiva ( 2005 ). iGomes et ah ( 2009 ). 



(2009 ). Als o, [Bhattj (12010) introduced the BS autoregressive 



Aslam et al.1 (120 10|) developed various acceptance sampling schemes 
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Figure 1: Plots of the density and hazard ratio function; (3 = 1. 



Little work has been published on the analysis of ce nsored data f or the BS distribution although 
censoring is common in reliability and survival studies. iRieckl (11995 ) consi dered parametric estima- 
tion for the BS distribution based on type-II censored samples, beng ( 2003 ) explored and compared 
different procedures to compute confidence interva ls for th e para meters and quanti les of the BS dis- 



tribution for complete and type-I censored data. iNg et al. 



(bood) 



and 



Wang et al.1 (120061) discussed 



the maximum likelihood estimation of the parameters of the BS distribution and proposed a simple 
bias-reduction method to reduce the bias of the maximum likelihoo d estimators ( MLEs ) based on 
type-II censored samples and random censoring, respectively. Also, iFrom and Li i2006l) discussed 
the estimation of the parameters of the BS distribution under censored samples. 

A particularly useful censoring mechanism is the type-II (right) censoring, which occurs when n 
items are placed on test and the experiment is terminated when the first m (m < n) items fail. Since 
life testing is often time consuming and expensive, type-II censoring may be used to reduce testing 
time and costs. This paper focus on testing inference for the BS distribution under type-II censored 
samples. 

A commonly used test procedure is the likelihoo d ratio test. A new criterion for testing hypothe- 



ses, referred to as the gradient test, was proposed by iTerrelll (|2002r) . Its statistic shares the same first 



order asymptotic properties with the likelihood ratio statistic. Like the likelihood ratio statistic and 
unlike the score and the Wald statistics, the gradient statistic does not require computation of the 
information matrix, neither observed nor expected, which is an advantage in problems involving cen- 
sored samples. The chief goal of our paper is to present Monte Carlo simulations in order to evaluate 
and compare the performance of the likelihood ratio and gradient tests for the BS distribution under 
type-II right censored samples. The other two tests will be briefly taken into account in our simulation 
study. To the best of our knowledge, there is no mention in the statistical literature on the use of the 



gradient test under censoring. 

The paper is organized as follows. Section [2] briefly describes the likelihood ratio and gradient 
tests. Section [3] presents inference on the parameters of the BS distribution under type-II right cen- 
sored samples. Numerical results are presented and discussed in SectionlU We evaluate and compare 
the performance of the likelihood ratio and gradient tests for testing hypotheses on the parameters of 
the BS distribution under type-II right censored samples. Applications are presented in Section [5J 
Finally, Section [6] closes the paper with some conclusions. 



2 Likelihood ratio and gradient tests 



Let £{6) denote the total log-likelihood function and consider the partition = (Oj , 0j) T , where the 
dimensions of 6i and 6 2 are q and p — q, respectively, i.e. 6 is a p- vector of unknown parameters. 
Let XJq denote the score function for 6. The partition for induces the corresponding partition: 
Uq = (Uj , Uj 2 ) T . Suppose the interest lies in testing the composite null hypothesis 



H : 2 







20 



against Hi : 2 ^ 9 2 o, where 2O is a specified vector. Hence, 0\ is a vector of nuisance parameters. 
The likelihood ratio (LR) and gradient (St) statistics for testing Ho versus Hi are given, respectively, 
by 

LR = 2{£(0) -£(0)}, S T = Uj(9-6), 



K^l > #20 



T \T 



denote the MLEs of = (6{ ,6jy under Hi and H 



with Ue 1 



where 6 = (Of ,0j) T and 6 = 
respectively, and U e = Uq(0) 
Uq 1 = 0, the gradient statistic can be written as St 
and St is Xp- q under Hq and Xp- q a> a noncentral chi-square distribution with p — q degrees of 
freedom and an appropriate noncentrality parameter A, under Hi. The null hypothesis is rejected for 



Uj 2 (0 2 



U 6l {0) and Ue 2 = U 02 (0). Since 
6 2 q). The limiting distribution of LR 



a given nominal level, 7 say, if the test statistic exceeds the upper 100(1 
distribution. 



7)% quantile of the % 



p-q 



Terrell 



(|2002[) points out that the gradient statistic "is not transparently non-negative, even though 



it must be so asymptotically." His Theorem 2 impl ies that if the log-likelihood function is concave 
and is differentiable at 6, then St > 0. Recently, 



Lemonte and Ferrari 



(2010a) derived the nonnull 



asymptotic expansion of the distribution of the gradient test statistic for a composite hypothesis under 
a sequence of local alternative hypotheses converging to the null hypothesis at rate n^ 1 ^ 2 . The authors 
also compared the local powers of the gradient, likelihood ratio, Wald a nd score tests. In general, no 
test is uniformly better than the others as far as local power is concerned. iLemonte and Ferraril (|20 10b) 
obtained the nonnull asymptotic expansions of the likelihood ratio, Wald, score and gradient statistics 
in Birnbaum-Saunders regression models when no censoring is present. An interesting finding is that, 
up to an error of order n^ 1 , the four tests have the same local power in this class of models. Their 
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simulation study evidenced that the score and the gradient tests perform better than the likelihood 
ratio and Wald tests in small and moderate-sized samples. 



3 Inference under type-II censored samples 

Let {ti, . . . ,t m } be an ordered type-II right censored random sample obtained from n units placed 
on a life-testing experiment wherein each unit has its lifetime following the BS distribution, with 
the largest (n — m) lifetimes having been censored. Let 6 = (a, 0) T be the parameter vector. For 
notational convenience, let 



H(z) 



1 - <£>(z) 



h(z) 



T(Z 



1/2 > 



a 



r = ± 
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where </>(•) is the standard normal density fun ction 



The likelihood function can be written as (|Balakrishnan and Cohenl . Il99l|) 

T{t*)\\n\{l-${v m )} n - 



L{0) = \ H^a,f3)ti 3/2 (U + f3)exp 



2a 2 



- , , m\(n — m)\ 
,i=i \ / j \ j 

where v m = p(t* m )/a. Thus, the log-likelihood function, except for a constant term, is given by 

m . m 

£{0) = £(a, fi) = m log{*(a, 0)} + £ logfe + 0) - — 2 £ r(t*) 



i=i 



i=i 



+ (n - m) log{l - $(v m )}. 



By taking partial derivatives of the log-likelihood function with respect to a and we obtain the 
components of the score vector, Uq = (U a , Up) T : 



a \ cr / cr 
x ' i=i 



'n — m)v m H(v r , 



a 



-HL _u 1 1 nf ,*2\ , (J^z m)h{t* m )H{v ri 

20 + 2^t l + + 2a 2 0^ P ^ ) + 

1=1 1=1 



2/3 



Setting U a and Up equal to zero yields the MLE 6 = (a, 0) T of0 = (a, 0) T . These equations cannot 

lem numerically via ite rative 



1999 



Press et al. 



2007) with 



be solved analytically and statistical software ca n be used to solve ti 
methods. The BFGS method (see, for example, iNocedal and Wright! 
analytical derivatives has been used for maximizing the log-likelihood function £(0). 

We now consider hypothesis testing on the parameters a and based on the likelihood ratio and 
gradient statistics. The interest lies in testing the null hypotheses 



"Hqo ■ oc = a 



(o) 



H 01 : = {o \ 
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which are tested against H 10 : a ^ and H n : (3 ^ (3^\ respectively. Here and (3^ are 
positive known scalars. For testing "H o> the likelihood ratio and gradient statistics are given by 

LR (a) = 2{£(p,a)-£(a®,P)}, S T(a) = U a {a-a®), 

where a and (3 are the unrestricted MLEs obtained from the maximization of £{d) under the alternative 
hypotheses, f3 is the restricted MLE obtained from the maximization of £(6) under H m and U a = 

U a (a(°\ (3). For testing Hqx, we have 

LR m = 2{£(j3,a) - i(a,^)}, S m = U p - /3 (0) ). 

Here, a is the restricted MLE obtained from the maximization of £{6) under the hypothesis "H i 
and Up = Up(a, /3(°'). The limiting distribution of these statistics is xl under the respective null 
hypothesis. In both cases, the null hypothesis is rejected if the chosen test statistic exceeds the upper 
100(1 — 7)% qua ntile of the y ? distribution. 



Theorem 3 in 



Terrell 



(|2002|) points out an important feature of the gradient test. It suggests that we 
can, in general, improve the approximation of the distribution of the gradient statistic by a chi- square 
distr ibution under the null hypothesis by using a less biased estimator. On the other hand, Theorem 
2 in 



Terrell 



(|2002f) indicates that if the estimate used in the gradient statistic is not the maximum 



likelihood, the non-negativity of the gradient statist ic is not g uaranteed. 



A bias-corrected estimator for a was derived in 



a 





( m\ 




. 1 + 2 H 1 -n). 


[ n 





Ng et al.l (|2006|) and is given by 
1 



The authors showed through Monte Carlo simulations that a is less biased than the original MLE, 
a, of a. Therefore, based on this estimator, we can define the adjusted gradient statistic = 
U a (a — a^) for testing the null hypothesis Hoo : ot = The limiting distribution of S^,^ is 
also x\ under Hqq. It is expected that the gradient test that uses the statistic S^, a ^ have better size 
performance than the gradient test based on the original gradient statistic Sr(a)- A disadvantage of 
the adjusted gradient statistic S^O) ^ s ma * ^ can ta ^ e on negative values. However, in order to avoid 
negative values we shall redefine S^,, a \ as 

S. 



T(a) 



max{0, U a {a-a {0) )}. 



An adjusted gradient statistic for testing the null hypothesis "H i : P — P^ is not considered here 
because a bias-corrected estimator for f3 is not available. In the next section, we shall present an 
extensive Monte Carlo study in order to evaluate and compare the performance of the tests presented 
in this section. 

4 Numerical results 



In this section we shall present an extensive Monte Carlo simulation study in which we evaluate and 
compare the finite sample performance of the likelihood ratio and gradient tests for testing hypothe- 
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ses on the parameters of the BS distribution under type-II right censored samples. We set the degree 
of censoring (d.o.c.) at 0(10)50%, the sample size at n = 20 and 40, and the shape parameter at 
a = 0.1, 0.3, 0.5, 0.75 and 1.0. Without loss of generality, the scale parameter j3 was kept fixed at 
1.0. The nominal levels of the tests were 7 = 10% and 5%. The number of Monte Carlo replications 
was 10000. All the Monte Carlo sim ulation experiments were performed using the Ox matrix pro- 
gramming language iDoornik , 2006 ). Ox is freely distributed for academic purposes and available at 
http : //www. doornik. com[ 

Table \T\ presents the null rejection rates (entries are percentages) of the tests which are based on 
the statistics LR^, S T ( a ) and ^(a) f° r testm g me nui l hypothesis V, o '■ ct = aS°\ Our main findings 
are as follows. First, for complete data without censoring (d.o.c. = 0%), the tests based on Sr( a ) and 
-S't(q) ^ l ess s ^ ze distorted than the test that uses LR( a y, in fact, they produce null rejection rates 
that are very close to the nominal levels in all the cases considered. In addition, the adjusted gradient 
test is slightly superior than the original gradient test. For example, for n = 20, = 0.5 and 
7 = 10%, the null rejection rates are 1 1.48% (LR {a) ), 9.65% (S T{a) ) and 9.81% (£* ( A Second, the 
size distortion of all tests increases with the d.o.c, the likelihood ratio test and the original gradient 
test displaying null rejection rates that are, respectively, greater and smaller than the nominal level. 
The adjusted gradient test is the less size distorted in most of the cases. For instance, for n = 20, 
a (°) = 0.75 and 7 = 5%, the null rejection rates are 6.17% (LR {a) ), 4.28% (S T(a) ) and 4.77% (5* (a) ) 
for d.o.c. = 10%, and 7.22% (LR {a) ), 4.15% (S T{a) ) and 4.96% (S* {a) ) for d.o.c. = 30%. Note that 
all the tests become less size distorted as the sample size increases, as expected. 

The null rejection rates of the tests which are based on LRm\ and St(b) for testing the null hy- 
pothesis "H i : /? = 1 are presented in Tabled For complete data without censoring (d.o.c. = 0%), 
the gradient test is less size distorted than the likelihood test in all the cases. For example, for n = 20, 
a = 0.5 and 7 = 10%, the null rejection rates are 1 1.98% (LR W ) and 10.53% (S T (j3)). Additionally, 
the size distortion of the tests increases with the degree of censoring and, similarly to what occurs for 
testing on the parameter a, the likelihood ratio test presented a liberal behavior, while the gradient 
test was conservative in the majority of the cases. For instance, when n = 20, a = 0.5 and 7 = 10%, 
the null rejection rates are 12.28% (LR {/3) ) and 9.88% (S T{/3) ) for d.o.c. = 20% and 13.57% (L% } ) 
and 9.47% (S T ^) for d.o.c. = 50%. It should be noticed that the gradient test is less size distorted 
than the likelihood test in all the cases considered. Again, the tests become less size distorted as the 
sample size increases. 

Table |3] contains the nonnull rejection rates (powers) of the tests. We set a = 0.5, (3 = 1, 
n = 80, 120 and 150, and the degree of censoring at 0%, 20% and 40%. The rejection rates were 
obtained under the alternative hypotheses H 10 : a = 8\ and H n : (3 = 5 2 , for different values of 5i 
and 82. The test that uses LR^ presented smaller powers for testing hypotheses on the parameter a 
than the tests which are based on the statistics Sr( a ) and S^, a y For example, when n = 80, d.o.c. = 
20% and 6t = 0.60, the nonnull rejection rates are 58.92% (LR (a) ), 62.02% (S T(a) ) and 65.53% 
(^r(a))- Additionally, the tests which use LR^ and S T ^ for testing hypothesis on the parameter 
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Table 1: Null rejection rates (%) for "H o : ot = a^. 











n = 


20 










n = 


40 






d.o.c. 




LR { , 


=0 


St( 


a) 


q* 




LR { 


a) 


ST(a) 


q* 


a) 


(%) 




10% 


5% 


10% 


5% 


10% 


5% 


10% 


5% 


10% 


5% 


10% 


5% 





0.1 


11.53 


6.07 


9.70 


4.25 


9.83 


4.52 


10.76 


5.23 


9.67 


4.82 


9.81 


5.16 




0.3 


11.50 


6.10 


9.70 


4.22 


9.78 


4.48 


10.74 


5.22 


9.69 


4.83 


9.83 


5.15 




0.5 


11.48 


6.15 


9.65 


4.21 


9.81 


4.45 


10.75 


5.24 


9.71 


4.80 


9.85 


5.12 




0.75 


11.55 


6.24 


9.58 


4.22 


9.78 


4.41 


10.77 


5.29 


9.74 


4.79 


9.83 


5.09 




1.0 


11.73 


6.25 


9.62 


4.26 


9.78 


4.46 


10.84 


5.36 


9.67 


4.81 


9.79 


5.10 


10 


0.1 


11.85 


6.12 


9.57 


4.29 


9.84 


4.90 


10.79 


5.70 


9.63 


4.61 


9.87 


4.79 




0.3 


11.86 


6.12 


9.55 


4.29 


9.82 


4.88 


10.77 


5.70 


9.62 


4.60 


9.85 


4.75 




0.5 


11.85 


6.15 


9.59 


4.29 


9.81 


4.85 


10.78 


5.67 


9.63 


4.56 


9.80 


4.73 




0.75 


11.96 


6.17 


9.59 


4.28 


9.82 


4.77 


10.76 


5.66 


9.59 


4.51 


9.88 


4.74 




1.0 


12.09 


6.19 


9.58 


4.35 


9.77 


4.73 


10.90 


5.69 


9.67 


4.51 


9.89 


4.73 


20 


0.1 


12.26 


6.50 


9.80 


4.10 


10.01 


5.09 


10.71 


5.30 


9.23 


4.36 


9.38 


4.63 




0.3 


12.24 


6.52 


9.74 


4.12 


9.96 


5.06 


10.77 


5.26 


9.24 


4.36 


9.29 


4.63 




0.5 


12.27 


6.50 


9.79 


4.16 


9.94 


5.04 


10.78 


5.25 


9.29 


4.37 


9.19 


4.63 




0.75 


12.26 


6.55 


9.68 


4.20 


9.91 


5.00 


10.82 


5.36 


9.32 


4.27 


9.33 


4.55 




1.0 


12.27 


6.64 


9.72 


4.22 


9.81 


4.93 


10.85 


5.36 


9.26 


4.28 


9.21 


4.52 


30 


0.1 


12.98 


6.98 


9.76 


4.30 


10.30 


5.12 


11.10 


5.55 


9.40 


4.43 


9.72 


4.64 




0.3 


13.04 


6.96 


9.74 


4.26 


10.17 


5.10 


11.11 


5.52 


9.40 


4.33 


9.69 


4.67 




0.5 


13.04 


7.06 


9.82 


4.28 


10.14 


5.04 


10.98 


5.50 


9.38 


4.30 


9.66 


4.73 




0.75 


13.04 


7.22 


9.81 


4.15 


9.96 


4.96 


10.97 


5.58 


9.34 


4.24 


9.56 


4.63 




1.0 


13.08 


7.26 


9.76 


4.19 


9.76 


4.91 


10.90 


5.59 


9.24 


4.24 


9.27 


4.60 


40 


0.1 


13.30 


7.23 


9.42 


4.27 


9.93 


5.15 


11.78 


5.89 


9.96 


4.79 


9.93 


5.38 




0.3 


13.31 


7.25 


9.40 


4.21 


9.95 


5.18 


11.75 


5.89 


9.98 


4.81 


9.95 


5.36 




0.5 


13.29 


7.12 


9.42 


4.14 


9.92 


5.13 


11.74 


5.93 


10.00 


4.75 


9.88 


5.26 




0.75 


13.25 


7.26 


9.32 


4.05 


9.90 


5.04 


11.73 


5.87 


9.83 


4.78 


9.82 


5.24 




1.0 


13.52 


7.32 


9.20 


3.95 


9.73 


4.95 


11.73 


5.95 


9.73 


4.68 


9.74 


5.09 


50 


0.1 


14.35 


7.99 


9.30 


3.95 


9.78 


4.91 


11.42 


6.02 


9.38 


4.57 


9.61 


5.10 




0.3 


14.37 


7.94 


9.24 


3.88 


9.70 


4.84 


11.39 


6.05 


9.42 


4.61 


9.64 


5.07 




0.5 


14.18 


7.92 


9.20 


3.89 


9.65 


4.85 


11.39 


6.12 


9.33 


4.55 


9.57 


5.03 




0.75 


14.37 


8.01 


9.04 


3.81 


9.56 


4.71 


11.42 


6.17 


9.39 


4.45 


9.39 


5.02 




1.0 


14.60 


8.05 


8.67 


3.58 


9.23 


4.59 


11.51 


6.23 


9.12 


4.33 


9.22 


4.91 
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Table 2: Null rejection rates (%) for Hoi '■ P — 1. 



n = 20 n = 40 



Li?( /3 ) 5 T ( /8 ) ^-R(/3) 'S't(^) 



(%) 


a 


10% 


5% 


10% 


5% 


10% 


5% 


10% 


5% 





0.1 


12.06 


6.43 


10.74 


4.98 


11.00 


5.43 


10.38 


4.92 




0.3 


11.95 


6.44 


10.57 


5.03 


11.02 


5.44 


10.37 


4.94 




0.5 


11.98 


6.52 


10.53 


4.94 


11.00 


5.50 


10.40 


4.95 




0.75 


11.78 


6.49 


10.33 


4.89 


11.14 


5.56 


10.26 


5.08 




1.0 


11.74 


6.74 


10.21 


4.93 


11.23 


5.81 


10.33 


5.08 


10 


0.1 


12.00 


6.44 


10.69 


4.70 


10.82 


5.59 


10.17 


4.76 




0.3 


12.01 


6.40 


10.40 


4.70 


10.80 


5.51 


10.16 


4.82 




0.5 


11.96 


6.40 


10.22 


4.63 


10.82 


5.45 


10.18 


4.81 




0.75 


11.91 


6.27 


9.89 


4.51 


11.02 


5.51 


10.15 


4.83 




1.0 


11.87 


6.29 


9.78 


4.60 


11.09 


5.64 


10.05 


4.80 


20 


0.1 


12.24 


6.44 


10.08 


4.57 


11.13 


5.55 


10.07 


4.75 




0.3 


12.31 


6.41 


10.00 


4.68 


11.14 


5.52 


9.96 


4.65 




0.5 


12.28 


6.36 


9.88 


4.52 


11.00 


5.55 


9.96 


4.56 




0.75 


12.32 


6.39 


9.72 


4.59 


11.01 


5.59 


9.81 


4.58 




1.0 


12.32 


6.53 


9.46 


4.44 


11.16 


5.47 


9.80 


4.61 


30 


0.1 


12.87 


6.95 


10.53 


4.66 


11.17 


5.60 


9.97 


4.55 




0.3 


12.92 


6.97 


10.23 


4.50 


11.18 


5.62 


9.95 


4.63 




0.5 


12.99 


6.99 


9.98 


4.39 


11.11 


5.63 


9.81 


4.57 




0.75 


13.01 


7.08 


9.70 


4.31 


11.01 


5.59 


9.65 


4.47 




1.0 


13.12 


6.91 


9.41 


4.39 


11.07 


5.69 


9.36 


4.45 


40 


0.1 


13.14 


7.08 


10.11 


4.38 


11.32 


5.84 


9.90 


4.45 




0.3 


13.14 


7.05 


9.67 


4.38 


11.32 


5.85 


9.65 


4.40 




0.5 


13.12 


7.05 


9.40 


4.26 


11.31 


5.85 


9.59 


4.34 




0.75 


13.10 


7.03 


9.05 


4.29 


11.22 


5.90 


9.52 


4.41 




1.0 


13.14 


6.99 


8.80 


4.39 


11.17 


5.92 


9.34 


4.24 


50 


0.1 


13.49 


7.85 


9.90 


4.22 


11.70 


5.98 


10.08 


4.49 




0.3 


13.52 


7.82 


9.67 


4.06 


11.71 


5.98 


9.88 


4.42 




0.5 


13.57 


7.75 


9.47 


4.24 


11.71 


5.97 


9.70 


4.41 




0.75 


13.54 


7.73 


8.81 


4.06 


11.75 


6.02 


9.38 


4.34 




1.0 


13.58 


7.69 


8.19 


4.23 


11.74 


6.03 


9.02 


4.33 
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(3 have similar powers, the gradient test being slightly more powerful than the likelihood test. For 
instance, when n = 150, d.o.c. = 0% and S 2 = 113, the nonnull rejection rates are 92.02% (LR^) 
and 92.33% (Sr(/3))- Note that the powers of the all the tests decrease as the degree of censoring 
increases. We also note that the powers of the tests increase with n and also with Si and 82, as 
expected. 

A comment on the values of the shape paramet er chosen for th e simu lation study, namely a = 
0.1, 0.3, 0.5, 0.75 and 1.0, is in order. Following Dupuis and Mills! (11998b. in most practical a ppli- 
cations, < a < 1; see, for inst ance, the applications c onsid ered in 



Leiva et al 



Sanhueza et al 



_ (120081) and 

(120081) . As showed in iHOyland and Rausandl (119941 . p. 46^8), this follows from the 
derivation of the BS distribution. As a result, most of the simulation experimen ts for the BS dis- 
tributi o n have been conducte d for values of a between and 1; s e e, for example. iDupuis and Mills 
Jl998l).lwuandWong| hoO^t) . kg et ali (I2006I) . Iwang et all (bood) . 



Cvsneiros et al. 



Lemonte et al. 



(2007 



2008) and 



(12008k Likewise, we chose values of a in the unit interval to run our simulations and 
to analyse the simulation findings in detail. We now briefly comment on simulation results for larger 
values of a. We conducted Monte Carlo experiments for various values of a > 1. In general, for 
testing on a and /3, the likelihood ratio test is oversized and the gradient test becomes conservative. 
However, the best performing test is the gradient test. For example, when n = 20, d.o.c. = 10%, 
7 = 10% and a = 1.8, the null rejection rates of the tests were 12.44% {LR {a) ), 9.77% (S T (a)) 
and 9.83% for testing H m : a = 1.8, and 11.90% (LR {/3) ) and 9.09% (S TW ) for testing 

Uoi:0 = 1. 

A natural question at this point is why the Wald and the Rao score tests were not included in the 
previous simulation experiments. Recall that the Wald and the score statistics involve the Fisher infor- 
mation matrix, which cannot be obtained for the BS distribution under type-II censoring. A common 
practice in such a case is to use the observed information in the place of the expected information. We 
followed this approach and ran various simulation experiments including the Wald and the score tests. 
Our general conclusion is that these tests cannot be recommended for the following reasons. First, 
the Wald and score tests were markedly oversized. For example, for n = 20, 7 = 10% and a = 0.75, 
the null rejection rates of the Wald and score tests for testing the null hypothesis "Hoi f3 = I are, 
respectively, 13.79% and 15.61% for d.o.c. = 10%, 16.73% and 15.44% for d.o.c. = 30% and 19.02% 
and 13.70% for d.o.c. = 50%. Second, in our simulations the inverse of the observed information 
matrix frequently produced negative standard errors in censored samples. 

In summary, the best performing test for a is the adjusted gradient test, i.e. the one that uses 
the bias-corrected estimator of a. As far as hypothesis testing on (3 is concerned, the gradient test 
performs better than the likelihood ratio test and should be preferred. 
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Table 3: Nonnull rejection rates (%): a = 0.5, j3 = 1, 7 = 10% and different sample sizes. 





d.o.c. 






a = 61 




W11 :/3 = 


62 


n 


(%) 


61 


LR( n \ 


St(c<) 


T(a) 


5 2 






80 





0.50 


10.71 


10.11 


10.25 


1.00 


10.53 


10.20 






0.54 


23.90 


25.84 


27.97 


1.04 


19.06 


19.50 






0.58 


56.16 


58.99 


62.18 


1.10 


54.60 


55.66 






0.60 


71.59 


73.97 


76.45 


1.13 


73.29 


73.96 




20 


0.50 


10.65 


9.55 


9.68 


1.00 


10.55 


10.10 






0.54 


19.34 


20.93 


23.62 


1.04 


18.67 


19.27 






0.58 


44.89 


48.51 


52.02 


1.10 


52.70 


54.15 






0.60 


58.92 


62.02 


65.53 


1.13 


70.86 


72.06 




40 


0.50 


10.93 


9.92 


9.72 


1.00 


10.30 


9.51 






0.54 


15.74 


17.20 


19.83 


1.04 


17.34 


18.47 






0.58 


33.78 


37.30 


40.92 


1.10 


47.69 


50.16 






0.60 


44.65 


48.93 


52.86 


1.13 


64.54 


67.06 


120 





0.50 


10.53 


10.17 


10.18 


1.00 


10.62 


10.37 






0.54 


30.34 


32.59 


34.90 


1.04 


22.74 


23.25 






0.58 


72.73 


74.76 


76.52 


1.10 


69.69 


70.35 






0.60 


86.61 


87.73 


88.78 


1.13 


87.01 


87.47 




20 


0.50 


10.64 


10.02 


10.14 


1.00 


10.19 


9.80 






0.54 


24.58 


26.43 


28.95 


1.04 


22.78 


23.46 






0.58 


60.06 


62.86 


65.65 


1.10 


68.10 


69.23 






0.60 


75.07 


77.04 


79.19 


1.13 


84.79 


85.37 




40 


0.50 


10.65 


9.90 


9.81 


1.00 


10.21 


9.65 






0.54 


19.02 


21.07 


23.71 


1.04 


20.63 


22.17 






0.58 


45.24 


48.49 


52.05 


1.10 


60.84 


62.71 






0.60 


59.85 


62.66 


65.77 


1.13 


79.08 


80.55 


150 





0.50 


9.98 


9.72 


9.76 


1.00 


10.15 


10.00 






0.54 


35.53 


38.04 


40.24 


1.04 


26.34 


26.97 






0.58 


80.99 


82.36 


83.63 


1.10 


77.31 


77.87 






0.60 


92.42 


93.24 


93.82 


1.13 


92.02 


92.33 




20 


0.50 


10.13 


9.94 


9.99 


1.00 


10.39 


10.15 






0.54 


28.30 


30.38 


33.01 


1.04 


25.52 


26.39 






0.58 


68.24 


70.56 


73.17 


1.10 


75.83 


76.92 






0.60 


83.46 


84.88 


86.39 


1.13 


91.42 


91.90 




40 


0.50 


10.51 


10.06 


9.85 


1.00 


10.33 


10.02 






0.54 


21.85 


23.99 


26.53 


1.04 


22.49 


23.86 






0.58 


53.25 


56.23 


59.28 


1.10 


68.41 


70.32 






0.60 


68.45 


70.74 


73.44 


1.13 


86.73 


87.78 
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5 Applications 



In this section, we shall illustrate the use of the likelihood ratio and gradient tests in three real data 
sets. The first data set, collected by members of the Instrument Development Unit of the Physical 
Research Staff, Boeing Airplane Company, are the number of cycles-to-failure of 101 strips of 6061- 
T6 aluminum sheeting, cut parallel to the direction of rolling. Each had been subjected to periodic 
loading with a frequen cy of 18 cycle s per s econd, and maximum stress of 31 kpsi. The data were 

(|2007|) . Suppose we are interested in testing H o : a = 0.15 
125 against H\\ : f3 ^ 125. These hypotheses were considered 



taken from Table 8 of 



Lemonte et al. 



ag ainst "Hm : a ^ -15 an d Hoi ■ ft 



in 



Cvsneiros et al. 



(2008). Assuming different censoring proportions, the observed values of the 
different test statistics and the corresponding p- values are given in Table HI Note that for complete 
data without censoring (m = 101) and when m = 80, at the 5% significance level, the likelihood 
ratio test does not reject the null hypothesis Hoo unlike the gradient tests. For the other values of 
m, the likelihood ratio and the gradient tests lead to the same conclusion. In all cases, the gradient 
tests provide stronger evidence against the null hypothesis (smaller p- values than the likelihood test). 
For testing the null hypothesis Hoi, the same decision is reached by the likelihood ratio and gradient 
tests, i.e. the null hypothesis Hoi is rejected at any usual significance level. We also note that as the 
censoring proportions increase, the p- values decrease. 



Table 4: Test statistics (p-values between parentheses); aluminum fatigue data. 







Inference on a 




Inference on j3 


rn 


LR( a ) 


Sr(a) 


°T(a) 




St(/3) 


101 


3.5771(0.0586) 


3.9841 (0.0459) 


4.3171(0.0377) 


9.4279 (0.0021) 


9.2402(0.0024) 


95 


2.8573(0.0910) 


3.1598(0.0755) 


3.4821(0.0620) 


9.4250(0.0021) 


9.2582(0.0023) 


90 


3.0826(0.0791) 


3.4342(0.0639) 


3.7969(0.0513) 


9.5167 (0.0020) 


9.3800(0.0022) 


80 


3.8361 (0.0502) 


4.3641 (0.0367) 


4.8300(0.0280) 


9.8999(0.0017) 


9.8573(0.0017) 


70 


2.8684(0.0903) 


3.2360(0.0720) 


3.6615(0.0557) 


9.4412(0.0021) 


9.5939(0.0020) 


60 


4.5172(0.0336) 


5.3218(0.0211) 


5.9240(0.0149) 


10.8407(0.0010) 


11.3070(0.0008) 


50 


4.0608(0.0439) 


4.8212(0.0281) 


5.4073(0.0201) 


10.3798(0.0013) 


11.3522(0.0008) 


40 


8.8234(0.0030) 


11.6943(0.0006) 


12.7236(0.0004) 


15.3808(0.0001) 


18.0787(0.0000) 



As a second application, we consider the data provided by Mc Cooll (1 19741) o n the the lifetime, 



Cohen et al. 



in hours, of 10 sustainers of a certain type. The data, used by 
of the three-parameter Weibull distribution, can be found at Table 5 of 
wish to test H o '■ a = 0.21 against H™ : a ^ 0.21 and % m : ft 



( 198 4j) as an illust ration 
i2007k We 



Lemonte et al. 



These hypotheses were considered in 



Cvsneiros et al. 



180 against Hu : p ^ 180. 
(|2008|) . Under different censoring proportions, 



the observed values of the different test statistics and the corresponding p- values are given in Tabled 
Note that for complete data (without censoring) the likelihood ratio test does not reject the null hy- 
pothesis Hoo at any usual significance level, whereas the original and adjusted gradient tests reject the 
null hypothesis Hoo at the 10% and 5% significance levels, respectively. Except for the case where 
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m = 10 (no censoring), the likelihood and gradient tests lead to the same conclusion. For testing the 
null hypothesis Hoi, unless m = 7, the same decision is reached by the likelihood ratio and gradient 
tests, i.e. the null hypothesis is rejected at the 10% significance level for m = 8, 9 and 10, and it not 
rejected at any usual significance level for m = 6. 



Table 5: Test statistics (p- values between parentheses); sustainers data. 







Inference on a 




Inference on (3 


m 


LR( a ) 




c* 
T(a) 


LR (P) 


St(/3) 


10 


2.1646(0.1412) 


2.7944(0.0946) 


4.0043 (0.0454) 


2.9417(0.0863) 


2.7580(0.0968) 


9 


0.0770(0.7814) 


0.0728(0.7873) 


0.0000(1.0000) 


3.2449(0.0716) 


2.9248 (0.0872) 


8 


0.3307 (0.5653) 


0.2911(0.5895) 


0.0000(1.0000) 


3.1616(0.0754) 


2.8499(0.0914) 


7 


0.6732(0.4119) 


0.5514(0.4578) 


0.1472(0.7013) 


2.9510(0.0858) 


2.7036(0.1001) 


6 


0.8471(0.3574) 


0.6620(0.4159) 


0.2234(0.6365) 


2.6797(0.1016) 


2.5463(0.1106) 



Finally, our third example uses the data given by iGuptal (|1952|) . which are the survival times (in 
days) of the first seven of a sample of ten mice after being inoculated with a culture of tuberculosis. 
The observations are: 41, 44, 46, 54, 55, 58 and 60. Suppose we are interested in testing "H o : a = 0.1 
against Hio : a ^ 0.1 and "H 01 : f3 = 54 against Hu : f3 ^ 54. We have LR( a ) = 1-7607 (p-value 
= 0.1845), S T ( a) = 2.3152 (p-value = 0.1281), 5* (a) = 3.7346 (p-value = 0.0533), LR m = 1.3710 
(p-value = 0.2416) and St(p) = 1.2054 (p-value = 0.2722). Note that the likelihood ratio and original 
gradient tests do not reject the null hypothesis T-Loo at any usual significance level, whereas the adjusted 
gradient test rejects the null hypothesis Hoo at the 10% level. Additionally, the null hypothesis 'Hoi is 
not rejected at any usual significance level based on the statistics LR^ and St(p). 



6 Concluding remarks 

The BS distribution is becoming increasingly popular in lifetime analyses and reliability studies. 
In this paper, we dealt with the issue of performing hypothesis testing concerning the parameters 
of this distribution under type-II right censored samples. We considered the likelihood ratio and 
a recently proposed test, the gradient test. As evidenced by our simulation results, in general the 
likelihood ratio test is oversized (liberal), whereas the gradient test is undersized (conservative) and 
less size distorted than the likelihood ratio test. We also considered an adjusted gradient statistic for 
testing hypotheses on the parameter a whose null distribution is more accurately approximated by the 
limiting null distribution than that of the original gradient test statistic. Overall, our numerical results 
favor the gradient tests for testing hypotheses on the parameters a and 0. The results in this paper 
are encouraging and motivate further work. We plan to investigate the performance of the gradient 
test as compared to that of the likelihood ratio test in the Birnbaum-Saunders regression model under 
various censoring mechanisms. 
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